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Abstract 

Superpositions of paraxial Laguerre-Gauss laser beam modes to generate optical potentials based 
on the optical dipole force are investigated theoretically. Parabolic and other monomial potentials 
with even powers, in cylindrical and spherical symmetry, with large diameters, can be synthesized. 
This superposition approach promises to help with high quality atom cloud manipulation and 
imaging. 
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I. INTRODUCTION 



The fields of atom [l[ and molecular optics 



m 



have developed considerably over recent 



years. Atomic beams that behave like laser beams can be created |4|: they are considered 
for direct deposition used in microscopy j^, precision metrology ol, and 



of quantum statistical effects 



m 



or the studies 
9|. 



such as the Hanbury Brown-Twiss effect 
The quality of the preparation of atomic ensembles, particularly in the cases of ultra- 
cold gases, is excellent, frequently as good as fundamental limits allow 



10 



.11 1. The same 

rr 

is not true of atom-optical components, for example atom-beam focussing [JJ suffers from 
the lack of lenses with large numerical apertures, see 12| and references therein. Similarly, 
wide parabolic optical dipole potentials for focussing and defocussing jol of atomic clouds 



at tolerable laser beam power 



12 



131 ] are currently unavailable. 



Here, the use of the optical dipole force for the application of wide high-quality aberration- 
free pulsed or stationary 2D and 3D isotropic potentials, primarily for the manipulation of 
cold gas clouds, is considered. In the regime of laser light far detuned from the atomic tran- 
sition the resultin g po tentials are conservative and their strength is proportional to the laser 
light intensity [ll. Il4|. It is shown that superpositions of Laguerre- Gaussian beams using 
spatial light modulators allow us to tailor laser beams in a suitable fashion [l^ to create 
wide cylindrical and spherical optical potentials improving control for the manipulation of 
atomic or molecular clouds. Harmonic potentials are considered in greatest detail because of 
their importance for focussing 12] and defocussing [ol, [3] in atom optics and atomic lithog- 
raphy l| or for feeding into waveguides [3] of atomic chips [l6]. The superposition approach 
can be extended to non-harmonic potentials with cylindrical and spherical symmetry, see 
section IIIII below. Non-harmonic potentials should facilitate dynamical studies of trapped 
atom clou ds [l7 |. perturbations, the formation of caustics, and aspects of quantum-classical 
behaviour 18 1. 

For focussing and defocussing of atomic beams many approaches have been investigated, 
see references in 12|. For focus sing and defocussing of cold gas clouds in two and three 
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21 



22| setups have been experimentally imple- 



dimensions magnetic focussing 

mented [l^. When miniaturizing such approaches care has to be taken that unwanted 



interactions with bulk media do not disturb trapped gas 



The superposition approach 



introduced here is meant to complement these and make the great flexibility of optical dipole 



2 



force approaches more accessible. 

In section |TT] we briefly recap the physics of Laguerre-Gauss modes LGp,i and the optical 
dipole force. We find that superpositions of members of the family {LGp i,p = 0, . . . , Nmax} 
are required for the generation of parabolic potentials, other values of orbital angular mo- 
mentum / can yield purely monomial potentials of order 21. We determine the expansion 
coefficients needed to create the desired optical potentials with cylindrical symmetry in sec- 
tion mil An important consequence of the superposition approach and the main motivation 
for this work is the ability to very considerably save laser power, see subsection IIII CI Wide 
aberration-free optical dipole force potentials become realizable which currently are infeasi- 
ble because of laser beam dilution, this has been recognized as a road block for some time: 
"In spite of numerous impressive achievements, using laser light interaction to tailor atom 
beams demands a high quality of the transverse collimation and this technique is difficult 
to scale with accessible laser powers" [l3|. 



The use of Laguerre-Gauss beams is not only experimentally well established 2J] they also 



form an appropriate basis for our analysis: limitations of the superposition approach arise 



because of mode-dispersion due to Gouy's phase |12l . l25l | . This phase is incorporated in the 
definition of Laguerre-Gauss modes, see Eq. ([1]); its effects are considered in subsections IIII Dl 
and lIV7\l We finally show how crossing two modulated beams with cylindrical symmetry can 
be used to form a spherically symmetrical potential in section [IV] and conclude in section |Vl 



II. LAGUERRE-GAUSS BEAMS 

Laguerre-Gauss modes are monochromatic, paraxial beam solutions 
and are defined as 
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Here, the monochromatic plane wave factor exp[i{kLZ—uJit)] is omitted |25l.l26l.l27l.l28l.l29|. 
L are generalized Laguerre polynomials, p and / the integer valued nodal, and angular mo- 
mentum numbers respectively, r = (p, z) is the position vector with the transverse coor- 
dinate vector p = {x,y)] the Kronecker-delta function Sqi refiects the fact that the modes 



with zero angular momentum have to be normahzed differently to the other modes. The 
transverse coordinates also parameterize the orbital angular momentum phase ^ via the 
relationship e*^ 



X 



ly 



261 ] ■ The frequency of the monochromatic laser gives rise to 



its wavenumber = uj^/c = 2tx / where A^, is the laser light's wavelength. The wave 
front radii R{z) = {z^ + zj^j/z, the beam radii w{z) = Wqa/I + / z% with the focal beam 

) = axctaji^z/ zr) are all 



radius wq = ^/XlZrJtt, and the longitudinal Gouy-phase shi: 



ts 



29|. 



parameterized by the beams' Rayleigh lengths zji |26l . 1271 . |28 

Correctly chosen superpositions of modes using spatial light modulators 

p 

p=l 
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30| 



(2) 



allows us to create superpositions ^ which give rise to parabolic intensity distributions. 



Following reference 28|] the use of [y-) polarized modes in Equation ([2])) yields an electric 
field which is polarized in the y-direction with a small contribution in the ^-direction due 
to the tilt of wave fronts off the beam axis (x, y, z are the unit- vectors and 3? stands for 
real-part) 



Ep{r- 1) = 3?{[y + z ic ^]e*(^^^-^*) 



}• 



(3) 



In keeping with the paraxial approximation of not overly focussed beams we neglect the 
transverse derivative in eq. ([3]). The associated time-averaged light intensity distribution 
then has the form 



/p(r)=eo {Ep{r,tf)^'-^ujl\^>p{r)\' . (4) 



A. Normalization, Intensity Scaling and Gradient Reduction 

With the normalized modes of Eq. ([1]) and assuming that the sum of the coefficients 
^ |cpp in Eq. ([2]) is normalized to unity we use the cross-sectional beam power normalization 



/oo POO 2 1*00 POO 2 
/ dxdy\^p{x,y,z)\^ = j dx dy Ip{x,y, z) = -Ip = 1 . (5) 
-oo J — oo oo J — oo 
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We note that the intensity /(x, z) of beams of fixed total power reduces inversely 
proportionally to their width wq in one direction, that is, their field amplitudes scale with 
WQ^^'^in X and y. Furthermore the field gradients diminish with Wq^. This implies that 
the effective curvature of the integrated laser light intensity | |V^p27rr(ir, responsible 
for atomic focussing, scales with Wq"^. We face an unfavourable quartic scaling with the 
beam width if we attempt to expand a laser beam transversally in p in order to widen the 
effective potential without weakening its power. Additionally, as we will show below, pure 
modes have small useful areas to generate the desiredpotentials, the combination of these 



two factors makes a pure mode approach unfeasible [12|, ll3|, |3ll]. It forces us to employ the 
mode superpositions studied here. Below, we discuss two approaches for the compensation 
of gradient weakening: by power compensation in subsection IIII A| and by beam waist 
narrowing in subsection IIIIBI 
B. Optical Dipole Force 

We assume that the interaction between atoms and the laser light is well described by 
a two-level scheme (excited state e and ground state g) in rotating wave approximation 
with effective atomic line width F and resonance frequency u = Uf. — ujg. This leads to the 
expression /(r) r'^/{2Is) = Q{r)'^ for the Rabi-frequency f2 as a function of the ratio of the 
local laser intensity I{r) and the transition's saturation intensity Is = 7r/icF/(3A'^) [l^. I32I. 
With sufficiently weak laser intensity I and sufficiently large detuning 5^ = cu^ — of the 
laser frequency ujl from the atomic transition frequency u, the AC-Stark shift gives rise to 
a conservative optical dipole potential which, to first order in I /Is, has the form [ssl 

8 6^ Is 16 d^Is 

This potential is modified due to detrimental spontaneous emission noise and light fluc- 
tuations. These tend to increase with increasing laser intensity but can be decrease d by 



increased detuning or through the use of more complicated optical level schemes 33 1 



Further discussion of their influences is beyond the scope of this paper. 



C. Potentials of order 21 



We now consider cylindrical atom-potentials with purely monomial modulation in the 
transverse direction oc p^'. A Taylor-expansion in p shows that LG- modes with angular 
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orbital momentum / depend in leading order on and every other higher order (p'"'"^"',n a 
positive integer). The Taylor coefficients of different modes are linearly independent of each 
other. Combining them into suitable superpositions created from mode families {LGp^i^p = 
0, . . . , Nmax} allows US to retain the leading and remove all higher order terms up to and 
including that of order p'+^^maa. Almost purely monomial potentials of order 21 therefore 
arise from such superpositions. 



III. CYLINDRICAL POTENTIALS 



The coefficients for superpositions using the mode family LGp^i are straightforward to 
determine through Gauss elimination. The ffist six superpositions {\l/p, P = 1, . . . , 6} yield 
the following (normalized X^^^^ |cppP = 1) coefficient matrix 
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FIG. 1: Amplitude coefficients cpp and probabilities Cp^ of Laguerre-Gauss superpositions ^'p of 
up to Nmax = S*'^ order modes {P = Nmax + 1 = 1, 2, . . . , 6) according to the coefficient matrix ([7]). 

Mode-superpositions extend the "useful" linear part of the field profile yielding wider 
parabolic intensity proffies. Figures [2] and [3] demonstrate that the parabolic part in the focal 
intensity profile of a superpositions field grows with the number of modes used. 
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FIG. 2: Left: Transverse electric field profile, Ep{x, 0,0), and Right: transverse intensity pro- 
files Ip{x, 0,0) at focal cross-section of Laguerre-Gauss beams comprising superpositions ^p, 
(P = 1,...,6), according to coefficient matrix d?]) (x-axis in units of focal beam radius wq, total 
cross-sectional beam power normalized to unity, (eo<^|^/2 set to unity), Rayleigh lengths zji kept 
constant). The vertical bars mark in the p 



confirming harmonic oscillator-scaling 



ot on the right are located at positions 0.57-\/2P -|- I-wqi 



25( 1 of the superposition beams' widths. 



A. Increased Beam Powers Compensate for Potentials' Widening 



If we increase the total cross-sectional beam power Ip for wider beam profiles according 
to the ratios of the modes' transverse derivatives, Ip = li\dx'^i{x,0,0)/dx'^p{x,0,0)\'^, 
the weakened gradient is power-compensated for by increased laser power. This way all 
superpositions give rise to potentials with equal strength, see Fig. [3l the necessary beam 
power increase to achieve this compensation is sketched in the inset of Fig. [31 



B. Decreased Rayleigh-Lengths Compensate for Potentials' Widening 

Alternatively to the beam-power increases just discussed, we can keep the total beam 
power for all beams equal and shrink the higher-order superposition-beams' Rayleigh lengths 
through increased beam focussing in the p-direction. This Rayleigh length-matching also 
allows us to compensate for the gradient reduction observed in Fig. O The laser intensity 
profiles for Rayleigh- matched superpositions are displayed in Fig. HI the filled-in areas in this 
figure are limited by the points dp, where each intensity curve deviates from the enveloping 
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FIG. 3: Focal intensity profiles I{x, 0) of Laguerre-Gauss beam superpositions comprising up to 
IQth Qj-(jgj- modes (same units as in Fig. [2]). For the superposition modes depicted in the left panel 
the total beam power has been adjusted such that all profiles have the same curvature at the origin 
as the dotted line parabola. The inset shows the necessary relative power increase, Ip/Ii, as a 
function of superposition order P (solid red line); it scales approximately like ^ ■ (inset: dotted 
black line). 



parabola (dotted line) by 0.74 percent. They delineate the useful areas of the potentials. 
This quality-criterion is adopted from Gallatin and Gould's work 3l| which showed that 
beyond a deviation of 0.74% spherical aberrations distort the atomic point-spread function 
of an imaged atomic beam too severely; for more details see references 12| and 311 ]. 



C. Power Savings 

The filled-in areas in Fig. H] represent the laser power fraction contributing to the atom 
potential in each case. Higher-order superpositions clearly allow us to use the laser power 
much more efficiently. Most of the laser power is wasted in the wings if the superposition 
approach is not employed. Additionally to the quantification of the useful area of the 
potentials delineated by the deviation points dp (see Fig. |4] and Table [T]). This waste is 
meaningfully quantified through the determination of the fraction of power Sp the laser 
beam contributes to the 'useful' part of the potential profile. We define it as the ratio of the 
laser energy contributing to the area between the deviation points \p\ < dp, in terms of the 
total laser power, namely 
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FIG. 4: Focal intensity profiles l{x,0,0) of Laguerre-Gauss beams comprising up to 16^^ 
order modes and their 0.74%-deviation marks dp, which lie at relative positions dp/di = 
1.00, 3.36, 9.96 from the origin, compare Table |T] (same units as in Fig. In contrast to 
Fig. [3] all superpositions have the same total beam power /, but the Rayleigh lengths zr have been 
readjusted such that all higher-order superpositions match up with curvature of the first mode 
case ^'i, see text. 

TABLE I: Potential Parameters dp and £-p, compare FigU] 
P 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 

dp/di 1.00 3.36 4.98 6.07 6.86 7.46 7.93 8.31 8.63 8.90 9.14 9.34 9.52 9.68 9.83 9.96 
£p [%] 0.0027 0.35 1.7 3.7 6.1 8.5 11 13 15 17 19 21 23 24 26 27 
£.p/£i 1 127 614 1362 2221 3103 3967 4791 5570 6306 6991 7639 8244 8816 9352 9858 




^ /o /o dr /p(pcos(r), p sin(r), 0) 
^ j^oo'^y j^^dx Ip{x,y,Q) 

Table [T] and Fig. O summarize and quantify our findings: Table [T] allows us to compare 
values for a single-mode atom potential, for which £i = 0.0027%, with the superposition 
approach. For example, compared to mode = LGq^i the relative power savings in case of 
superposition \l/i6 is 9858, this translates into a power utilization of Siq = 0.0027% x 9858 = 
27%. In general the details of this behaviour depend on the chosen quality criterion but 
the underlying scaling is straightforward to derive. The useful fraction of the laser beam is 
proportional to a 2D integral over the intensity and therefore grows with the fourth power 
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of the position of the deviation mark EpjEi = {dp/di)^, for example Siq/Si = (diQ/di)^ ~ 
9.96^ ^ 9858. 

£p/£i 




5 10 15 20 25 30 

P 



FIG. 5: Relative power savings £p/£i as a function of superposition mode number P. 



D. Limitations due to Mode Dispersion 

Different order modes carry different Gouy-phase factors which leads to mode dispersion 
in the focal region 25j. The associated change in intensity distribution is illustrated in 
Fig. [6l An interesting case is the creation of spherically symmetrical potentials using two 
beams with cylindrical symmetry and equal foci orthogonally crossing each other. In this 
case the width of one beam dictates the area along the beam axis of the other beam that is 
used. This implies constraints on the degrees of permissible focussing which we investigate 
in the following section. 



IV. SPHERICAL POTENTIALS 



We now want to investigate the scenario for the generation of spherically symmetrical 
potentials. If an identical copy of the laser beam that travels along the z-axis is additionally 
sent along the x-axis such that their crossed confocal configuration leads to the simultaneous 
application of two cylindrical potentials a spherical potential is applied to the gas cloud. 
The laser beams have to be sufficiently detuned from each other in order to avoid harmful 
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FIG. 6: Intensity profile Iq{x, 0, z) of Laguerre-Gauss superposition from two perspectives. The 
mode dispersive effects of Gouy's pliase cliange tlie intensity distribution along the beam axis. Note 
that the x-axis is displayed in terms of the focal beam radius wq whereas the z-axis is represented 
in terms of the Rayleigh length zr. 



interference despite their spatial overlap [3J]. They moreover have to be elliptically stretched 
in the y-direction, by a factor \/2 in the parabolic case or by 2^/(20 in the case of a monomial 
potential oc p^', otherwise adding up their intensities would remove the desired isotropy of 
the spherical potential. 

A. Limitations due to Mode Dispersion 

Gouy's phase (f){z) = aTctan{z / zn) ^ z/zr, varies strongest near the beam focus and 



introduces relative phases between the modes within each beam [25|. If the beam is very 
strongly focussed (small value of zr) the dephasing away from the focus 2; = is so rapid 
that non-linear aberrations degrade the desired linear field profile. In other words, a lower 
limit for the Rayleigh lengths Zminl-P) as a function of the number of used modes P has to 
be determined in order to guarantee moderate dephasing. Whereas the absolute values for 
this lower limit are hard to derive from first principles, we can still work out the correct 
scaling with the maximal mode number P: 

The electric field is proportional to the superposition of the modes including the Gouy- 
phase factors; this can be approximated by Ep oc J2p=o ^p+i-^^p^^*^"^ ~ J2p=o Cp+i-^G'p(l + 
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FIG. 7: The left panel illustrates the behavior of the relative deviation of the intensity distribution 
AI from zero as it approaches the 0.74%-deviation marks (top and bottom grid). Here, A/30 is 
shown for the crossed configuration of two laser beams travelling along z and x-axis respectively. 
The value of the Rayleigh length at which we find that the oscillatory behaviour of AI along 
a constant radial perimeter just exhausts the upper and lower limits set by the deviation marks 
allows us to determine the associated value of Zmin- The latter is plotted as a function of maximum 
mode number, in the middle panel (the filled in blue area is the forbidden area of too tightly 
focussed beams). The values of z^[^(P) in turn determine the position of the turning points 
0.57-V2PTT-W0 (top green line, compare Fig. [2]), the position of the deviation-points dp (middle 
blue line), and the focal beam radii iL'o(-Zmm(/')) (lower red line), depicted in terms of the laser's 
wavelength in the right panel. 

2ipz/zji). The expansion coefficients are positive and the wave functions are real at the 
focus z = 0. Since the first order term is purely imaginary the intensity has to depend on z 
quadratically: Ip{x,y,z) = lp{x,y,0) ■ [1 + ^Dp + 0{z^)]. The deviation term Dp has a 
complicated dependence on the number of modes, but, containing the square of sums of the 
form X]^=o^ pCp+iLGp, is roughly proportional to P^. When we consider the relative deviation 
of the intensity profile near the focus from the focal intensity distribution, AI = -^^^j^^y^, 
we find AJp oc ^ ■ P^. Additionally, we know that the widths of the superpositions scale 



For constant relative intensity deviations AJp this implies const. = ■ P^ or zr oc P^^'^. 
A numerical investigation, see Fig. H confirms z^i^{P) = 2.2 ■ ■ P^/^ as a good estimate 
for a lower bound on zp. This relationship has been checked numerically and holds for 
7 < P < 30. There is no reason to believe deviations might occur for values of P > 30, but 



roughly like those of the harmonic oscillator 25|, see Fig. |2l namely z oc V2P + 1 ^ V2VP. 
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for small values of P the assumptions used in the derivation of the scaling law do not hold 
accurately, see Fig. [2j Instead, the expression ZminP = 15 ■ Al ■ P^/^ gives a better estimate 
for ZrainP the range of < P < 6. These lower limits for zr imply that the beam focus 
is several wavelengths wide and a posteriori confirms that the paraxial approximations hold 
for all cases discussed here, since the largest beam opening angle conforming with the lower 
limits presented here turns out to be roughly 7° for superposition \l/2- 



V. CONCLUSIONS 



For a possible experimental implementation of focussing fields repulsive (blue-detuned) 
optical potentials with a dark center are probably most suited since they minimize detri- 
mental spontaneous emission noise. For an expansion field analogously red detuned optical 
potentials should be used. If fine-tuning is considered one will probably also have to 
revisit the approximations underlying Eq. (jl]) and Eq. (jSD; such considerations are beyond 
the scope of this paper. 

The techniques for the coherent superposition of laser modes have been experimentally 



demonstrated, see e.g. references 2J, |30|| and citations therein. We have found here that 



using the mode-superposition approach allows for very considerable laser power savings and 
potentials can be made wider than is possible with pure modes. We come to the conclusion 
that for the design of atomic potentials, based on the optical dipole force, it is possible and 
necessary to coherently superpose suitable laser modes in order to create wide high quality 
parabolic potentials. 
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